Spin freezing and dynamics in CasCo^-zMn^Oe (x ~ 0.95) investigated with implanted 
muons: disorder in the anisotropic next-nearest neighbor Ising model 
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We present a muon-spin relaxation investigation of the Ising chain magnet Ca3Co2- : rMn a; 06 
(x ~ 0.95). We find dynamic spin fluctuations persisting down to the lowest measured temperature 
of 1.6 K. The previously observed transition at around T « 18 K is interpreted as a subtle change 
in dynamics for a minority of the spins coupling to the muon that we interpret as spins locking 
into clusters. The dynamics of this fraction of spins freeze below a temperature Tsf ~ 8 K, while a 
majority of spins continue to fluctuate. An explanation of the low temperature behavior is suggested 
in terms of the predictions of the anisotropic next-nearest-neighbor Ising model. 
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The magnetic chain multiferroic CasC^-zMn^Og 
(a; « 1) has been the subject of considerable re- 
cent investigational^. This material is based on 
the Ising spin chain magnet CasC^Og, with (close 
to) half of the cobalt ions replaced with manganese^. 
The observation of up- up-down-down (TTIi) order in 
this system has led to the proposal that, at low 
temperatures, Ca3Co2-;cMn x 06 may be described by 
the anisotropic next-nearest-neighbor Ising (ANNNI) 
model 1 . This modet 8 ^ ' 10 : 11 describes Ising spins on 
a three-dimensional lattice in which, along one direc- 
tion, there is nearest-neighbor ferromagnetic exchange 
(Jfm) and next-nearest-neighbor antiferromagnetic ex- 
change (Jafm)- For |Jafm/<7fm| > 1/2 the ground 
state magnetic order is of the TTIi type. As temper- 
ature is increased from T = the magnetic behav- 
ior is determined by the existence of domain wall soli- 
tons which separate regions with different commensurate 
AFM spin arrangement o 9 ' 10 . Although a continuum de- 
scription of the ANNNI model predicts an infinity of high 
order commensurate AFM phases (known as the devil's 
staircase) a description in terms of a discrete Hamilto- 
nian shows the possibility of metastable states of ran- 
domly pinned solitons. In magnetic systems, these so- 
called "chaotic states" are expected to lead to frozen-in 
disorder or spin glass- like behavior 8 . Here we present 
an investigation of the low temperature static and dy- 
namic magnetism in Ca3Co2_ x Mn x 06 (x ks 1) that we 
have observed at a local level using muon-spin relaxation 
(/i + SR). We find that the low temperature magnetic state 
of Ca3Co2- a; Mn 2 ;06 is reached through a complex freez- 
ing out of dynamic processes and we conjecture that the 
existence of chaotic states provides an explanation for the 
disordered magnetism and persistent dynamics that we 
observe at low temperature. 

Ca3Co2- K Mn a; 06 is formed from chains of magnetic 
ions arranged along the c-axis in alternating oxygen cages 
of face-shared trigonal prisms and octahedra. Mn 4+ ions 
preferentially occupy the octahedral sites while the trig- 



onal prisms are occupied by Co 2+ ions. The magnetic 
chains form a triangular lattice in the ab plane separated 
by Ca 2+ ions. While it is agreed that 3d 3 Mn 4+ is in the 
5 = 3/2 high-spin configuration, the spin state of 3d 7 
Co 2+ has been questioned. Although fits to magnetic 
neutron diffraction data 1 - suggest a low spin S = 1/2 
state, electronic structure calculations 3 ^ and x-ray ab- 
sorption spectroscopic 3 - favor the high spin S = '3/2 
state. The neutron diffraction and magnetic suscepti- 
bility measurements 1 - indicate that below T-q w 18 K 
spins align along the c-axis, adopting the TTII config- 
uration with Mn and Co ordered moments of 0.66/xb and 
1.93/iB respectively. The considerable width of the mag- 
netic Bragg peaks suggests that this is not a state of 
true long-range order (LRO), but rather represents the 
locking- in of spins into finite sized domains. Taken with 
the Ising-like character of the magnetic ions, the obser- 
vation of TTII order is suggestive that this material can 
be described as a realization of the ANNNI model, or 
its extension— to the case of chains of unequal Ising 
spins. It is also notable that the inversion symmetry 
breaking of TTII order, along with the alternating charge 
order, results in magnetism driven ferroelectricity in this 
material^. 

/x + SR has proven useful in elucidating the static and 
dynamic properties of Ising systems^ 3 - including the par- 
ent compound CasC^Oe^ 1 ^- In a /^ + SR experiment, 
spin-polarized positive muons are stopped in a target 
sample, where the muon usually occupies an interstitial 
position in the crystal. The observed property in the ex- 
periment is the time evolution of the muon spin polariza- 
tion, the behavior of which depends on the local magnetic 
field at the muon site, and which is proportional to the 
positron asymmetry function A(t). Measurements were 
made at the Swiss Muon Source (S/xS) using the GPS in- 
strument. A polycrystalline sample of Ca3Co2-zMn. E 06 
with x = 0.95 (similar composition to Ref. 1), was packed 
in a Ag foil packet and mounted on a Ag plate in a helium 
cryostat. 
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FIG. 1: ZF /^+SR spectra measured at (a) T = 110 K, (b) 14.6 K and (c) 4.13 K. Solid lines are fits to Eq. JTJ. Inset to (b): the 
parameter y from Eq. ([1} shows evidence of static freezing around T ~ 8 K. Inset to (c): Spectra measured at 1.6 K, showing 
a small, heavily damped oscillation. The solid line is a guide to the eye. 



Spectra measured for Ca3C01.05Mno.95Og are shown in 
Fig. [0 For temperatures T > 10 K [Fig. Ufa) and (b)] 
we observe purely relaxing asymmetry spectra. We do 
not observe oscillations in the asymmetry, which usually 
signal the presence of long-range magnetic order (LRO). 
For T < 8 K the form of the spectra alters, most no- 
tably with an increase in the non-relaxing baseline of 
the asymmetry. There is also a small undulation in the 
asymmetry that emerges at early times in high statis- 
tics spectra measured at the lowest temperatures [inset 
FigHJc)]- This feature, which may represent a low am- 
plitude, heavily damped oscillation displayed little tem- 
perature dependence and was too indistinct to be fitted 
systematically. Instead, the spectra are best modelled 
over the entire temperature range by a sum of asymme- 
try functions, one with a large relaxation rate Ai and the 
other with a smaller rate A2: 

A(t) = A hg + Ao(pe~ x ^ (1) 

+{1 _ rt[ | + (l _| )e -«]), 

where Aq represents the amplitude of the signal arising 
from the sample and the term A\, g represents a temper- 
ature independent, nonrelaxing background from those 
muons that stop in the sample holder or cryostat tails. 
The parameter y — above ~ 10 K but takes nonzero 
values at low temperatures (see below). The amplitude 
p was found to be p = 0.75 across the entire measured 
temperature regime. Exponential relaxation is often ex- 
pected in cases where dynamic fluctuations in the local 
magnetic field at the muon site represent the dominant 
relaxation process^ and it is likely that this mechanism 
is at work in Ca3Co2_ a; Mn 2: 06. Fast field fluctuations 
lead to relaxation rates that vary as A; oc 7^(-B?)r, where 
7^ is the muon gyromagnetic ratio, Bi is the local mag- 
netic field at the ith muon site and r is a fluctuation 



rate. In the present case, the interpretation of dynamic 
fluctuations is supported by the observation that applied 
longitudinal magnetic fields of up to 0.6 T do not decou- 
ple the relaxation, as would be expected for relaxation 
from static field inhomogeneities^. 

The coexistence of two relaxation rates, Ai and A2, 
implies the existence of two classes of spatially separate 
muon sites in Ca3Co2- x Mn x 06. In general, these sites 
might differ in the width of the field distribution (B 2 ), or 
the fluctuation time at each site may be different. (It is 
not possible to have a single class of site with two corre- 
lation times giving rise to two different relaxation rates; 
the shorter time will always dominate, giving the smaller 
relaxation rate^.) The first class of muon site (with am- 
plitude p and relaxation rate Ai) accounts for 75% of 
the muon sites, with the second, with amplitude (1 — p), 
accounting for the remaining w 25% of muon sites. As 
there is no evidence (from our measurments or previous 
worki) for phase separation in this system, it is proba- 
ble that both classes arise from the intrinsic behavior of 
the bulk of the material. Possibilities for this include the 
coupling of each class of site preferentially to one or other 
of the magnetic cations (i.e. one class sensitive to fields 
arising from Co 2+ and the other sensitive to Mn 4+ ), as is 
the case in of X3V2O8 (X=Ni,Co) 2 °, or that the classes 
of muon site are sensitive to different components of the 
same spin distribution as in GeN^O^^L. In the latter 
case, a system with local site anisotropy might, for ex- 
ample, give rise to relaxation timescales for longitudinal 
and transverse fluctuations that could be quite different. 
If different muons sites were selectively sensitive to lon- 
gitudinal or transverse components, dynamics with two 
timescales may arise from the same spin site. 

The temperature dependence of Ai and A2 resulting 
from fits of the data to Eq. JT]) are shown in Fig. [5] 
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A| at = 88(3) MHz, U x 
and A| at 



Both relaxation rates show the same trend of behavior 
with their magnitude increasing as the temperature is re- 
duced, followed by saturation of the relaxation rate at a 
constant value below ss 30 K. This trend is sometimes 
seen in the /x + SR of complex systems with low tem- 
perature dynamics 22,2 ^ and its origin is not completely 
understood. It can arise due to the existence of two 
relaxation channels, one of which is strongly T depen- 
dent with a correlation time r s while the other shows 
little variation with T and has a correlation time r w . 
As noted above, in the presence of two competing relax- 
ation processes, that with the shorter correlation time 
wins out, giving the smaller relaxation rate. At high 
temperature, therefore, we have t s (T) <C r w which re- 
sults in a strongly T dependent relaxation which we can 
fit phenomenologically by Ai = Ci exp([/i/T). At low 
temperatures, where r w <C t s (T), we have A^ ~ Af at , re- 
sulting in a phenomenological fitting function l/Ai(T) = 
1/Af at + l/[Ciexp(Ui/T)]. Fitting this function to our 
data allows us to parameterize the relaxation rates with 
275(5) K for the fast relaxation 
2.2(2) MHz, U 2 ~ 270(20) K for the slowly 
relaxing component. Given the similarity between U\ 
and U2, it is likely that similar relaxation processes are 
at work at both sites, sharing similar correlation times. 
This would imply that the widths of the magnetic field 
distributions differ for the two classes of site and are in 
the ratio ((B 2 ) / (£f)) 1/2 « (AfVAf*) 1 / 2 « 6. We note 
further that, with the large activation energies involved, 
it is unlikely that the temperature dependent behavior 
is related to dynamics predicted by the ANNNI model. 
This is not surprising, since models of this sort provide a 
low energy description of real magnetic systems. Rather, 
it is likely that the observed behavior arises due to the 
single ion anisotropies of the magnetic ions. For free 
Co 2+ , for example, the spin orbit coupling constant^ 
is —270 K, which resembles our energy scale U. This 
is suggestive that the energy scale for the temperature 
dependent contribution to the relaxation rates is set by 
fluctuations between spin-orbit split components of the 
low-lying electronic states of the magnetic Co 2+ ions. It 
should be expected that Mn 4+ will affect the magnetic 
field distribution at the muon sites making the true sit- 
uation more complex. 

In addition to this general trend reflecting the dynam- 
ics of the system, the smaller relaxing component, with 
relaxation rate A2 , shows additional features at low tem- 
perature (although it is not obvious why these features 
are not seen in Ai). For the material in this regime, 
where the ANNNI model provides a good effective de- 
scription of the magnetic behavior, the spin fluctuations 
will include excitations involving the domain wall soli- 
tons described by the ANNNI model^, along with diffu- 
sive modes resulting from coupling of the Ising spins to 
non magnetic degrees of freedom. At around T B ~ 18 K 
we see a small peak in the relaxation rate A2 which occurs 
near the temperature below which (broadened) magnetic 
Bragg peaks appear in neutron diffraction^. The absence 
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FIG. 2: (a) Evolution of the relaxation rates Ai (open circles) 
and A2 (closed circles) as a function of temperature in ZF. 
(b) Detail from (a) showing departure of A2 from the trend 
around the transition at Tb and the spin freezing transition 
Tsf- Inset: Temperature evolution of the frequency of the 
small oscillatory feature seen below 8 K at early times. 



of a change at Tb in the form of the signal or its ampli- 
tudes indicates that this is not a transition to LRO or to 
a static magnetic state. Rather, it is more likely to repre- 
sent the freezing out of one or more relaxation processes, 
likely to involve the free motion of domain walls, leading 
to the formation of poorly correlated clusters of ordered 
spins. This is consistent with the observation of broad- 
ened peaks in the neutron diffraction 1 signalling that the 
order that gives rise to the magnetic Bragg peaks is not 
truly long range. 

The more significant change in the slowly relaxing com- 
ponent of the muon asymmetry occurs below 8 K. Here 
the parameter y in Eq. ([T]) increases on cooling from 
y = to y = 1 at the lowest measured temperatures, 
as shown in Fig. [T] indicating a transition to a static 
local field distribution at low T. Such a distribution, 
whether ordered or disordered, will only dephase those 
muon-spin components that lie perpendicular to the ini- 
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tial muon-spin polarization direction, expected to be 2/3 
of the total polarization for a powder sample, leaving the 
other 1/3 polarized. Although the magnetic field dis- 
tribution experienced by these muons below T m 8 K 
is static on the muon timescale, it is unlikely that the 
system locks into true LRO. Instead of the oscillations 
that would usually be observed in the presence of LRO, 
we see only a very small minimum in the asymmetry at 
early times [inset Fig.[lfc)] in the minority component of 
our spectra measured below 8 K. This may originate from 
Kubo-Toyabe-like relaxation 1 - 7 - typical of a static ensem- 
ble of disordered moments or may be a highly damped 
oscillation due to static order occurring over only a short 
length scale. Assuming the latter interpretation, the in- 
set to Fig. HJb) shows the frequency v{— j fl Bi/2ir) of this 
feature^. The static disordered magnetism observed in 
our experiments, taken together with the peak measured 
in the x" susceptibility and broad maximum in the heat 
capacity observed around 8 Ki, points to a freezing of 
domain walls to form spin clusters, resulting in a static 
disordered or glassy magnetic system— . There are also, 
of course, still significant slow dynamic fluctuations in 
the system, measured by the majority asymmetry com- 
ponent with amplitude p. 

The observed disorder may be explained in terms of the 
predictions of the ANNNI model^. Although the devil's 
staircase of commensurate spin structures is predicted 
from a continuum treatment of the model, a treatment 
based on a discrete Hamiltonian predicts the existence of 
an array of "chaotic" states at a higher energy than the 



devil's staircase, comprising a random array of domain 
wall solitons and antisolitons. Such states are metastable, 
but are separated from the true Devil's staircase of stable 
states by relatively high energy barriers, making it impos- 
sible for a system to relax into the devil's staircase states 
in a finite amount of time 8 -. Instead of LRO, the chaotic 
state possesses an intrinsic randomness which means it 
may be described as spin glass-like. This picture there- 
fore provides an explanation within the framework of the 
ANNNI model for the disordered magnetic state that we 
observe in CaC01.05Mno.95O6. We note, however, that 
additionals factors not considered in this simple model 
may contribute to the glassy character, including the pos- 
sibility of clustering of the Mn 4+ and Co 2+ ions and the 
effect of complex interchain exchange interactions 27 . 

Finally we note that it is possible that the behavior ob- 
served in CaC01.05Mno.95O6 reflects the existence of both 
series spin relaxation processes (where the freezing of one 
relaxation process allows another to freeze at a lower 
temperature) and also coexistent parallel relaxation pro- 
cesses which persist independently. Both of these types of 
process have been advanced to explain dynamics in glassy 
systems 2 - 8 - and it is likely that they are important in a 
complex dynamic system such as Ca3C01.05Mno.95O6- 
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